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Constraints on R n gravity from precession of orbits of S2-like stars 
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We study some possible observational signatures of R" gravity at Galactic scales and how these 
signatures could be used for constraining this type of f(R) gravity. For that purpose, we performed 
two-body simulations in R n gravity potential and analyzed the obtained trajectories of S2-like stars 
around Galactic center, as well as resulting parameter space of R n gravity potential. Here, we 
discuss the constraints on the R" gravity which can be obtained from the observations of orbits 
of S2-like stars with the present and next generations of large telescopes. We make comparison 
between the theoretical results and observations. Our results show that the most probable value 
for the parameter r c in R n gravity potential in the case of S2-like stars is ~100 AU, while the 
universal parameter ft is close to 0.01. Also, the R n gravity potential induces the precession of S2- 
like stars orbit in opposite direction with respect to General Relativity, therefore, such a behavior of 
orbits qualitatively is similar to a behavior of Newtonian orbits with a bulk distribution of matter 
(including a stellar cluster and dark matter distributions). 

PACS numbers: 04.50.Kd, 04.80.Cc, 04.25.Nx, 04.50.-h 



I. INTRODUCTION 



Power-law fourth-order theories of gravity have been 
proposed like alternative approaches to Newtonian grav- 
ity [3, Hi- In this paper we study possible application 
of R n gravity on Galactic scales, for explaining observed 
precession of orbits of S-stars, as well as weather these 
observations could be used for constraining this type of 
f(R) gravity 0. 

S-stars are the bright stars which move around the 
massive black hole in the center of our Galaxy • For 
one of them, called S2, there are some observational indi- 
cations that this orbit deviates from the Keplerian case 
due to relativistic procession Q. Besides, an extended 
dark mass which probably exists in the Galactic center, 
could also contribute to pericenter precessing of the S2, 
but in the opposite direction 0, @] ■ Progress in monitor- 
ing bright stars near the Galactic Center have been made 
recently Q. With the Keck 10 m telescope, the several 
stars orbiting the black hole in Galactic Center have been 
monitored, and in some cases almost entire orbits, as, for 
example, that of the S2 star, have been observed, allow- 
ing an unprecedented description of the Galactic Center 
region The astrometric limit for S2 star orbit is to- 
day around 10 mas and within that limit one can not say 
for sure that S2 star orbit really deviates from the New- 
tonian case. In the future, it will be possible to measure 
the positions of some stars with astrometric errors sev- 
eral times smaller than errors of current observations and 
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that is why we will consider here even smaller astrometric 
limits. 

Capozziello et al. [2[ investigated the possibility that 
the observed flatness of the rotation curves of spiral 
galaxies is not evidence for the existence of dark mat- 
ter (DM) haloes, but rather a signal of the breakdown 
of General Relativity (GR). They found a very good 
agreement between the theoretical rotation curves and 
the data using only stellar disc and interstellar gas when 
the slope n of the gravity Lagrangian is set to the value n 
— 3.5 (giving fi — 0.817), obtained by fitting the Type la 
supernova Hubble diagram with the assumed power-law 
f(R) model and without dark matter Q. 

Frigerio Martins and Salucci Q have also investigated 
the possibility of fitting the rotation curves of spiral 
galaxies with the power-law fourth-order theory of grav- 
ity, without the need for dark matter. They show that, in 
general, the power law f(R) version could fit the obser- 
vations well, with reasonable values for the mass model. 

Recently, gravitational microlensing has been investi- 
gated in the framework of the weak field limit of fourth 
order gravity theory Q . The solar system data (i.e. plan- 
etary periods) and light bending due to microlensing can 
be used to put strong constraints on the parameters of 
this class of gravity theories. In paper Q it was found 
that these parameters must be very close to those corre- 
sponding to the Newtonian limit of the theory. In paper 
[1C| | the authors discuss the constraints that can be ob- 
tained from the orbit analysis of stars (as S2 and S16) 
moving inside the DM concentration. In particular, con- 
sideration of the S2 star apoastron shift may allow im- 
proving limits on the DM mass and size. 

Rubilar and Eckart [ll[ investigated the properties of 
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stellar orbits close to central mass and the corresponding 
connection with current and (near) future observational 
capabilities. They showed that the orbital precession can 
occur due to relativistic effects, resulting in a prograde 
shift, and due to a possible extended mass distribution, 
producing a retrograde shift. Both, prograde relativistic 
and retrograde Newtonian periastron shifts will result in 
rosette shaped orbits. Weinberg et al. [l2[ discuss phys- 
ical experiments achievable via the monitoring of stellar 
dynamics near the massive black hole at the Galactic 
Center with a diffraction-limited, next-generation, ex- 
tremely large telescope (ELT). They demonstrate that 
the lowest order relativistic effects, such as the prograde 
precession, will be detectable if the astrometric precision 
become less then 0.5 mas. 

In this paper we continue to investigate constraints on 
the parameters of this class of gravity theories using S2- 
like star orbits under uncertainty of 10 mas. In Section 
§2 the type of used gravitational potential is given. In 
Section §3 we present the S2-like stars orbits, gravity pa- 
rameters and angles of orbital precession, and also com- 
pared theoretical results with observations. The main 
conclusions are pointed out in §4. 



II. THEORY 

R" gravity belongs to power-law fourth-order theories 
of gravity obtained by replacing the scalar curvature R 
with f(R) = foR n in the gravity Lagrangian [l|, [2j. As 
a result, in the weak field limit (l3j . the gravitational 
potential is found to be [H, : 



$ (r) = - 



GM 
~ 1 h r 



, T 
1 + I — 



(1) 



where r c is an arbitrary parameter, depending on the 
typical scale of the considered system and /3 is a universal 
parameter: 



12n 2 - In - 1 - V36n 4 + 12n 3 - 83n 2 + 50n + 1 
6n 2 - 4n + 2 ' 

This formula corresponds to a modification of the grav- 
ity action in the form: 



A = 



^g{f (R) + L m ), 



(3) 



where f(R) is a generic function of the Ricci scalar cur- 
vature and L m is the standard matter Lagrangian. 

For n — 1 and (3 — the R n potential reduces to the 
Newtonian one, as expected. Parameter j3 controls the 
shape of the correction term and is related to n which is 
part of the gravity Lagrangian. Since it is the same for 
all gravitating systems, as a consequence, (3 must be the 



same for all of them and therefore it is universal param- 
eter 0. The parameter r c is the scalelength parameter 
and is related to the boundary conditions and the mass 
of the system Q • 



III. RESULTS 

A. Orbits of S2-like stars and parameters of R n 
gravity 

In order to study the effects of R n gravity on motion 
of the star S2, we performed two-body calculations of its 
orbit in the R n potential (Eq.([T])) during two periods. 
We assumed the following input parameters taken from 
the paper of Zakharov et al. [Kj: orbital eccentricity 
of S2-like star e = 0.87, major semi-axis a — 919 AU, 
mass of S2-like star M* = 1 M Q , mass of central black 
hole Mbh = 3.4 xl0 6 M Q (where M Q is solar mass) and 
orbital period of S2-like star is 15 years. We calculated 
S2-like star orbit during two periods using Newtonian 
and R n potentials. We also investigated the constraints 
on the parameters (3 and r c for which the deviations be- 
tween the S2-like stars orbits in the R n gravity poten- 
tial (Eq.([T])) and its Keplerian orbit will stay within the 
maximum precision of the current instruments (about 10 
mas), during one orbital period. 

In Fig. [T]we presented trajectory of S2-like star around 
massive black hole in R n gravity (blue solid line) and in 
Newtonian gravity (red dashed line) for r c = 100 AU and 
for the following nine values of parameter (3: 0.005, 0.01, 
0.015, 0.02, 0.025, 0.03, 0.035, 0.04, 0.0475. The black 
hole is assumed to be located at coordinate origin. We 
fixed a value of parameter r c on 100 AU because this 
value corresponds to maximal value of parameter /3 in 
the parameter space (see Fig. [3]) and varied values of 
parameter j3. All 9 presented orbits fulfill request that 
R n orbit and corresponding Newtonian orbit differ less 
then 10 mas, (i.e. within the maximum precision of the 
current observations) during one orbital period. We can 
see that if parameter f3 increases R n orbit differs more 
from the corresponding Newtonian orbit since the preces- 
sion angle becomes larger. This indicates that the value 
of f3 should be small, as inferred from Solar system data 
@, and in contrast to the value (3 = 0.817 (obtained 
by Q which gives excellent agreement between theoret- 
ical and observed rotation curves). In the future, with 
improvements in observational facilities the precision on 
constrains on values of parameters /3 and r c will increase, 
as well as the accuracy of the S2 orbit. 

The corresponding distances between the S2-like star 
and black hole as a function of time for the same values 
of parameters r c and (3 as in the Fig. [T]are presented in 
Fig. [2 There is an additional requirement on parameter 
space: period of S2-like star orbit has to remain «15 ± 
0.2 yr. Like in previous case, with increasing observa- 
tional accuracy of period the precision on constraints on 
values of parameters /3 and r c will also increase. 
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FIG. 1: The orbits of an S2-like star around massive black hole in R n gravity (blue solid line) and in Newtonian gravity (red 
dashed line) for r c = 100 AU and for these nine values of parameter /3: 0.005, 0.01, 0.015, 0.02, 0.025, 0.03, 0.035, 0.04, 0.0475. 
The black hole with mass Mbh = 3.4 x 10 6 M© is assumed to be located at coordinate origin, and mass of S2-like star is taken 
to be M* = 1 M e . 



In Fig. |3] we presented the parameter space for R n 
gravity under constrain that, during one orbital period, 
S2-like star orbits under R n gravity differ less than e from 
their orbits under Newtonian gravity for 10 values of pa- 
rameter e: 0.001, 0.002, 0.003, 0.004, 0.005, 0.006, 0.007, 
0.008, 0.009 and 0".01. For e = CT.01, we can see that 
the maximal value of (3 is 0.0475, and the corresponding 
r c is 100 AU. That is why we investigate combinations 
P < 0.0475 and r c = 100 AU. This study is important 
because with improvements in observational facilities the 
precision of e will increase. 

In Fig. [4] we presented the parameter space for R n 
gravity under constraint that, during one orbital period, 
S2-like star orbits in R n gravity differ less than e from the 
corresponding orbits in Newtonian gravity for 12 values 



of parameter e: 0.01, 0.02, 0.03, 0.04, 0.05, 0.06, 0.07, 
0.08, 0.09, 0.1, 0.11 and 0"T2. Although these values are 
unrealistic, since they are larger then the current obser- 
vational resolutions, they can be used to obtain tendency 
of P and r c dependences on s. 

We want to find how the constraints on values of pa- 
rameters P and r c will change with increasing of e. From 
Figs. [3] and H] we can see that smaller values of e will sig- 
nificantly reduce parameter space of R n gravity within 
the given precision. 

In Fig. [SJleft and right) we presented dependence of 
the maximal value of parameter P versus e for < e < 
0".3 and < e < 0".l, respectively. We can see that 
in the region of interest the dependence of the maximal 
value of parameter P versus values of e is almost strictly 
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FIG. 2: The distances between the S2-like star and black hole as a function of time for the same values of parameters r c and 
P as in the Fig. Q] 



linear. 

In Fig. [Bfleft and right) we presented dependence of 
the value of parameter r™ ax , that correspond to maximal 
value of parameter /3 (e.g. f3 max ), versus e for < e < 
0".3 and < e < 0".l, respectively. With decrease of e 
the value of parameter r™ ax increase and in the region 
of interest the value of r™ ax is near 100 AU. 

The trajectories of S2-like star around massive black 
hole in R n gravity (blue solid line) and in Newtonian 
gravity (red dashed line) are presented in Fig. [7] (left and 
right) for r c = 100 AU and = 0.02 during 0.8 and 10 
periods, respectively. We can see that the precession of 
S2-like star orbit is in the clockwise direction in the case 
when the revolution of S2-like star is in counter clockwise 
direction. 

In Figure[5]we show calculated S2-like star orbits for 11 



periods, assuming r c — 100 AU and these nine values of 
parameter /3: 0.005, 0.01, 0.015, 0.02, 0.025, 0.03, 0.035, 
0.04, 0.0475. We adopted value r c = 100 AU because it 
corresponds to the largest allowed range of parameter /3. 
We obtained that the R n gravity causes periastron shifts 
which result in rosette shaped orbits. 



B. Angle of orbital precession in R n gravity 

Talmadge et al. [IH used post-Newtonian formalism of 
metric theories of gravity in order to calculate perihelion 
precession in a potential which deviates from Newtonian 
potential only slightly. Adkins and McDonnell [l5[ cal- 
culated the precession of Keplerian orbits under the in- 
fluence of arbitrary central force perturbations. For some 
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FIG. 3: The parameter space for R n gravity under the constraint that, during one orbital period, S2-like star orbits 
gravity differ less than e from the corresponding orbits in Newtonian gravity, for the following 10 values of parameter e: 
0.002, 0.003, 0.004, 0.005, 0.006, 0.007, 0.008, 0.009 and 0".01. 
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examples including the Yukawa potential they presented 
the results using hypergeometric functions. Schmidt [l6j 
calculated the perihelion precession of nearly circular or- 
bits in a central potential which has a form of modified 
Newtonian potential. Since the S2-like star orbit is very 
eccentric, and hence non-circular, we used approach from 
[lH to obtain analytical expression for precession angle 
in R n gravity. 

Assuming a potential which does not differ significantly 
from Newtonian potential, in our case R n gravitational 
potential, we will derive formula for precession angle of 
the modified orbit, during one orbital period. First step 



is to derive perturbing potential from: 

V(r) = $ (r) - <f N (r) ; ( r ) = . (4) 

r 

Obtained perturbing potential is of the form: 




and it can be used for calculating the precession angle 
according to the equation (30) from paper [15| : 
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(6) 



where r is related to z via: r 



1 



By differenti- 



ating the perturbing potential V{z) and substituting its 
derivative and expression for the semilatus rectum of the 
orbital ellipse (L — a (l — e 2 )) in above equation ©, we 
obtain: 



X2Fl ,l±i,^ ;2;e > 



(7) 



where a-F i is hypergeometric function. The graphical pre- 
sentation of the precession angle A8 for S2-like star orbit 
as a function of (3 is given in Fig. [9] (black solid line). 
From this figure it can be seen that A9 is negative for all 
values of /3 between and 1, which are of interest in the 
case of S2-like star orbit. 
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FIG. 5: The dependence of the maximal value of parameter fj on precision e ranging from to 0".3 (left) and from to 0".l 
(right). 
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FIG. 7: The orbits of S2-like star around massive black hole in R n gravity (blue solid line) and in Newtonian gravity (red 
dashed line) for r c = 100 AU and /3 = 0.02 during 0.8 periods (left) and 10 periods (right). 



Exact expression ([7]) is inappropriate for practical ap- in Taylor's series over f3, up to the first order, leads to 
plications. However, it can be easily approximated for the following expression for precession angle: 
/?«0 and (3 w 1. In case of (3 m expansion of Eq. ([7]) 
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FIG. 8: The orbital precession of an S2-like star around massive black hole located at coordinate origin in R" gravity for r c 
100 AU and these nine values of parameter /3: 0.005, 0.01, 0.015, 0.02, 0.025, 0.03, 0.035, 0.04, 0.0475. 



A0 = 



TT rad /3 (VT~^ - 1) 



(8) 



180°^ ( % /T _ ? - 1) 



Above expression in the case of S2-like star orbit is 
presented in Fig. [9] as a blue dash-dotted line. Similarly, 
expansion of Eq. (J7J in power series for f3 sa 1, leads to 
the following expression for precession angle (red dotted 
line in Fig. [H]): 



A0 



TT rad a (f3 - 1) _ i + e 2) 



180°a {ft - 1) (VT^ - 1 + e 2 ) 



(9) 



One can expect that in general precession angle de- 
pends on semimajor axis and eccentricity of the orbit 
(see e.g. Iorio & Ruggiero 2008 [17]), as well as on both 
potential parameters j3 and r c . It is indeed case for /3 ~ 
1 in Eq. (0). But as it can be seen from formula ©, 
the precession angle in the case when j3 is small (j3 ps 0) 
depends only on eccentricity and universal constant /3 
itself. 

In order to test if the approximation from Eq. ((5J) is 




FIG. 9: The exact expression for precession angle A9 in R n 
gravity (in degrees) as a function of the parameter /3 (black 
solid line) and two its approximations: for /3 ~ (blue dash- 
dotted line) and for /3 ~ 1 (red dotted line). The other pa- 
rameters correspond to the case of S2-like star: a = 919 AU, 
e = 0.87 and r c = 100 AU. 

satisfactory in case of S2-like star, we derived its preces- 
sion angle in two ways: 

• analytically from the approximative formula ([5]) 

• numerically from calculated orbits presented in Fig. 

m 

Comparison of the obtained precession angles by these 
two methods is presented in Table HI As it can be seen 
from this table, the approximative formula can be 
used for estimating the precession angle for all values of 
P from Fig. E 

Above analysis indicates that R n gravity results with 
retrograde shift of S2-like star orbit. Rubilar and Eckart 
PH showed that the orbital precession can be due to 
relativistic effects, resulting in a prograde shift, or due 
to a extended mass distribution, producing a retrograde 
shift. We can conclude that perturbing potential V(r) 
has a similar effect as extended mass distribution, since 
it produces a retrograde orbital shift. 

Since the precession has negative direction, as in the 
case of extended mass distribution, the obtained results 
are useful for testing if the precession due to extended 
dark matter enclosed into orbit of S2-like star could be 
also explained by R n gravity. If this is possible, it will ex- 
clude the need for dark matter hypothesis. Therefore, if 
future and more precise observations of bright stars near 
the Galactic Center will show a precession in the nega- 
tive direction, we have to conclude that the phenomenon 
could be caused by bulk distributions of stellar cluster 
or/and dark matter in classical Newtonian (GR) gravity 
or by R n gravity. On the other hand, if there is no devia- 
tions from Newtonian (GR) trajectories with an accuracy 
of observations, one could put constraints on stellar clus- 
ter and dark matter distributions and on parameters of 
R n gravity if we adopt the theory to fit observational 
data. 



9 



TABLE I: The numerically calculated by computer simulation 
(A6°) and analytically calculated from Eq. (O (A# c ) values of 
precession angle of S2-like star orbit (in degrees) as a function 
of universal constant /3 of R n gravity with parameter r c — 100 
AU. 



p 


A6»° 


A8 C 


0.005 


-0.602 


-0.604 


0.01 


-1.203 


-1.209 


0.015 


-1.802 


-1.816 


0.02 


-2.400 


-2.425 


0.025 


-2.997 


-3.035 


0.03 


-3.592 


-3.647 


0.035 


-4.186 


-4.261 


0.040 


-4.779 


-4.876 


0.045 


-5.666 


-5.493 



C. Comparison between the theoretical results and 
observations 

Here we compare the obtained theoretical results for 
S2-like star orbits in the R n potential with two inde- 
pendent sets of observations of the S2 star, obtained 
by New Technology Telescope/ Very Large Telescope 
(NTT/VLT), as well as by Keck telescope [see Fig. 1 
in|J], which are publicly available as the supplementary 
on-line data to the electronic version of the paper [J. 
However, all the above two-body simulations in R n grav- 
ity potential resulted with the true orbits of S2-like stars, 
i.e. the simulated positions of S2-like stars presented in 
Figs. [TJ [3J [7] and |S] are in their orbital planes. There- 
fore, in order to compare them with observed positions, 
the first step is to project them to the observer's sky 
plane, i.e. to calculate the corresponding apparent or- 
bits. From the theory of binary stars it is well known 
that any point (a;, y) on the true orbit could be projected 
into the point (x c , y c ) on the apparent orbit according to 
[see e.g.Hil: 

X° = l\X + hy, y c = m\x + m-iy, (10) 

where the expressions for li,l2,mi and m-i depend on 
three orbital elements (£1 - longitude of the ascending 
node, oj - longitude of pericenter and i - inclination) 
[ijj [2(| • One should take into account that in the case of 
orbital precession w is a function of time, and therefore 
should be in general treated as an adjustable parameter 
during the fitting procedure. However, the previously 
mentioned theoretical and observational results showed 
that in the case of S2-like stars this precession is most 
likely very small, and hence we assumed w as a con- 
stant when projecting true positions to their correspond- 
ing apparent values. For that purpose we used the fol- 
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FIG. 10: The fitted orbit of S2 star around massive black hole in R n gravity for r c — 100 AU and j3 = 0.01 (black solid 
lines in both panels). The NTT/VLT astrometric observations are presented in the left panel by blue circles, while the Keck 
measurements are denoted by red circles in the right panel. 
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FIG. 11: Comparison between the fitted (black solid lines) and measured (open circles) distances of the S2 star from black hole 
in the case of NNT/VLT (left) and Keck (right) observations. 



lowing Keplerian orbital elements from [4(: i = 134°. 87, 
f2 = 226°. 53 and uj = 64°. 98. Besides, our previous theo- 
retical results indicated that the most likely value of the 
scale parameter is r c ss 100 AU, and we adopted that 
value in order to reduce the number of free parameters 
when fitting the observations. 

We fitted the observed orbits of S2 star using the fol- 
lowing procedure: 

1. initial values for S2 star true position (xo,yo)i or- 
bital velocity (£o,2/o) an0 - ^ ne parameter (3 of R n 
gravity potential are specified; 



2. 



the positions (xi,yi) of the S2 star along its true 
orbit are calculated at the observed epochs using 
two-body simulations in the R n gravity potential, 
assuming that distance to the S2 star is = 8.3 
kpc and mass of central black hole Mbh = 4.4 
xl0 6 M Q 0; 

the corresponding positions (x°,yf) along the ap- 
parent orbit are calculated using the expression 

JED; 



3. 



4. the root mean square (O — C) goodness of fit is es- 
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FIG. 12: Comparison between the fitted (black solid line) and 
measured radial velocities for the S2 star. Measured velocities 
are labeled with open circles (VLT data) and open rhombuses 
(Keck data). 



timated according the following expression: 



(O - C), 



\ 



N 

E 

i=l 



(s? - x\f + (y? - ytr 



2N 



where (x°, y°) is the i-th observed position, (xf, yf) 
is the corresponding calculated position, and N is 
the number of observations; 

5. the values of the input parameters are varied and 
the procedure is repeated until the minimum of 



(O - C) r 



is reached. 



The best fit is obtained for the following small value 
of the universal constant: f3 = 0.01, in which case the 
corresponding precession is around —1° (see Table |IJ. In 
Fig. [TUlwe present two comparisons between the obtained 
best fit orbit for /3 = 0.01 in the R n gravity potential and 
the positions of S2 star observed by NTT/ VLT (left) and 
Keck (right). The corresponding calculated distances of 
S2 star from massive black hole are shown in Fig. QTJ As- 
trometric data for the S2 star orbit are presented by blue 
dots (NTT/ VLT measurements) and by red dots (Keck 
measurements). As one can see from these figures, there 
is a good agreement between the theoretical orbit and 
NTT/VLT observations. In case of Keck measurements, 
we had to move the origin of the coordinate system with 
respect to the both axes for 5 mas, in order to get reason- 
able fit. We made this correction following the suggestion 
from where it was necessary in order to combine the 
two data sets. After that we also achieved the satisfying 
agreement between the same fitted orbit and the both 
NTT/VLT and Keck data sets, in spite the fact that 
both groups obtained slightly different orbital elements, 
distance to the S2 star, as well as mass of central black 

hole 0,1 m. 

In order to obtain the orbital elements of S2 star both, 
NTT/VLT and Keck groups, fitted their observations 



with Keplerian orbits, but at the same time, they had 
to allow that the position of the center of mass varies as 
a function of time, i.e. they had to introduce the orbital 
precession. As it can be seen from Fig. [TU1 the orbit of 
S2 is not closed in vicinity of its apocenter, which clearly 
shows that the orbital precession is a natural consequence 
of R n gravity. Moreover, by comparing the arcs of orbit 
near the apocenter with the corresponding results pre- 
sented in Fig. 1 from [J], one can see that their curvatures 
are different, which indicates the opposite directions of 
precession in these two cases. Therefore, the future more 
precise observations of S2 star positions near its apocen- 
ter could have a decisive role in verifying or disproving 
the validity of R n gravity near the Galactic Center. 

We also made a comparison between the fitted and 
measured radial velocities for the S2 star (see Fig. [T2l . 
The well known expression for radial velocity in polar 
coordinates r and 9 is [see e.g. Il9j: 



Vrad 



sin i sin 



• + lo) ■ r + r cos(6» + u) -9 



(11) 



However, we used the corresponding expression in rectan- 
gular coordinates x = r cos 9 and y — r sin 9 to calculate 
the fitted radial velocities: 



Vrad — 



\Jx 2 + y 2 



[sm(9 + u>) ■ (xx + yy)+ 



(12) 



+ cos(# + uj) ■ (xy — yx)] 



y 

where 9 = arctan— . As it can be seen from Fig. Q21 the 

x 

agreement between our theoretical predictions and the 
observations is also satisfactory. 

Although the both observational sets indicate that the 
orbit of S2 star most likely is not a Keplerian one, the 
nowadays astrometric limit of around 10 mas is not suf- 
ficient to unambiguously confirm such claim. We hope 
that in the future, it will be possible to measure the stel- 
lar positions with astrometric errors several times smaller 
than errors of current observations. 



IV. CONCLUSIONS 

In this paper S2-like star orbit has been investigated 
in the framework of fourth order gravity theory. Using 
the observed positions of S2 star we put new constraints 
on the parameters of this class of gravity theories. We 
confirmed that these parameters must be very close to 
those corresponding to the Newtonian limit of the the- 
ory. For parameter /3 approaching to zero, we recover 
the value of the Keplerian orbit for S2 star. Also, we 
performed two-body calculations of its orbit in the R n 
potential. The obtained results showed that, in contrast 
to General Relativity, R n gravity gives retrograde direc- 
tion of the precession of the S2 orbit, like in the case 
when it is caused by an extended matter concentration 
in Newtonian potential. 
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Despite the excellent agreement between theoretical 
and observed rotation curves obtained by Capozzicllo 
and coworkers Q for R n parameter /3 (the slope n of 
the gravity Lagrangian is set to the value n — 3.5 giv- 
ing (3 — 0.817), our findings indicate that for e = 0".01 
maximal value of j3 is 0.0475, i.e. j3 is less or equal than 
0.0475, and our fitting indicated that optimal value for f3 
is around 0.01. Therefore, R n gravity in this form may 
not represent a good candidate to solve both the dark en- 
ergy problem on cosmological scales and the dark matter 
one on galactic scales using the same value of parame- 
ter /?. But this theory has its own benefits in explaining 
orbits of the stars and solar system data. 

For today astrometric limit of around 10 mas for S2 
star orbit, within that limit one can not say for sure that 
S2 star orbit really deviates from the Newtonian case, i.e. 
we have to stress that at the moment observations are 
in agreement with the Newtonian point-like potential for 
the Galactic Center. Therefore the observations and their 
theoretical analysis give us one of the best cases to discuss 
departures from the standard GR plus stellar cluster and 
dark matter (as it was done in our papers and papers 
of other authors) or to analyze an opportunity to get 



constraints on alternative theories observing trajectories 
of S2-like stars. The newest astrometric data for the star 
S2 of NTT/VLT measurements and Keck measurements 
show the Keplerian orbit fits for the respective data set, 
do not yield closed ellipses. Maybe this represents really 
small deviation of S2 star orbit from the Newtonian case 
and for more sure conclusion we need astrometric errors 
several times smaller than these errors, but we compared 
these data with S2 star orbit obtained using R n gravity 
potential. 

We can conclude that additional term in R n grav- 
ity compared to Newtonian gravity has a similar effect 
like extended mass distribution and produce a retrograde 
shift, that results in rosette shaped orbits. 
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